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Future changes in the number of dry days per year can either reinforce or counteract projected increases in 
daily precipitation intensity as the climate warms. We analyze climate model projected changes in the 
number of dry days using 28 coupled global climate models from the Coupled Model Intercomparison 
Project, version 5 (CMIP5). We find that the Mediterranean Sea region, parts of Central and South America, 
and western Indonesia could experience up to 30 more dry days per year by the end of this century. We 
illustrate how changes in the number of dry days and the precipitation intensity on precipitating days 
combine to produce changes in annual precipitation, and show that over much of the subtropics the change 
in number of dry days dominates the annual changes in precipitation and accounts for a large part of the 
change in interannual precipitation variability. 

Future drying in subtropical subsidence regions such as Mediterranean Europe is a recognized consequence 
in projections of anthropogenically forced climate change, but previous analyses of future precipitation have 
focused primarily on changes at annual or monthly time scales 13 . However, people and ecosystems experi- 
ence weather at the daily time scale so it is important to also anticipate changes in daily weather. Additionally, the 
frequency of dry days affects regional hydrology and ecosystems and potentially influences tourism, agriculture, 
and transportation. Previously it has not been possible to examine projected daily changes in precipitation in a 
large number of models 1,3 , but the CMIP5 archive contains daily precipitation fields from nearly 30 models, 
making such an analysis possible. By analyzing precipitation regimes at the daily scale we can now show the key 
role that trends in the occurrence of dry days play in regional changes in annual precipitation. 

Global projections of annual mean precipitation include robust responses to greenhouse gas (GHG) forcing in 
many areas, especially over tropical, subtropical and high latitude regions 3 5 . This response includes increases 
over historically wet areas of the deep tropics and high latitudes and decreases over semi-arid regions in the 
subtropics. The latter is a particular concern for society, as water resources are often scarce in these areas 
already 6 8 . Moreover, fewer precipitating days in dry areas should lead to increased interannual variability in 
annual total precipitation, exacerbating the erratic hydro-climates of arid, semi-arid and Mediterranean-type 
climate regimes. 

Our goals in this work are to answer such questions as how model simulations of the occurrence of dry days 
compare to observations and what the model-projected changes in dry day occurrence by the end of the century 
are. We demonstrate how changing numbers of days that are dry and changing precipitation intensity on days 
that are wet together determine annual precipitation changes, and show the effects of increasing numbers of dry 
days in some areas has on interannual precipitation variability. 

It is useful to place the current work on dry days in the context of previous work examining the physical 
mechanisms of future precipitation change. Thermodynamical mechanisms include the consequences of the 
Clausius-Clapeyron equation-governed increase in atmospheric water vapor concentration and transport in a 
warmer climate 4,9 . Dynamical mechanisms are related to poleward shifts of storm tracks and related cyclone/ 
anticyclone trajectories, baroclinic zones, and jets as descending Hadley cell branches and subtropical dry zones 
expand poleward and as midlatitude westerlies adjust to a reduced equator-to-pole temperature gradient 8,1012 . In 
some regions, only one of these mechanisms dominates; for example, Ref. 13 indicates that projected precipitation 
declines found between subtropical minima and midlatitude precipitation maxima arise from the poleward 
expansion of dry zones rather than from thermodynamical mechanisms. In other areas, the mechanisms com- 
pound or compete to yield projected overall precipitation changes. Mediterranean climates, for example, receive 
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Figure 1 | Dry day frequency (days/year) from (a) the CMIP5 multi-model ensemble average over the period of 1960-1989; (b) GPCP Daily 
observations averaged over the period 1997-2012. Maps were produced using NCAR Command Language (NCL; Ref. 30). 



the majority of their infrequent precipitation from the equatorward 
margin of midlatitude cyclones. Even a slight poleward shift in the 
storm-track can strongly decrease the probability of daily precipita- 
tion. At the same time, well-positioned atmospheric rivers are pro- 
jected to tap into and episodically deliver more moisture to some 
subtropical locales in a warmer climate, for example in California 14 . 
This work does not focus on these shifting mechanisms per se, but 
rather on the net effect climate change has on the occurrence of dry 
day and on its consequences for regional hydro-climates. 

Future increases in precipitation intensity are likely to be an 
important aspect of climate change, since warming will tend to 
accelerate the overall hydrological cycle, intensifying the wet 
extremes 5,15,16 in addition to the dry extremes in some areas. 
Several studies have quantified the impact of climate change on the 
wet tails of high frequency precipitation distributions 15,1719 , but stud- 
ies of the dry tails of regional precipitation distributions have up to 
now mostly focused on long-term drought, dry spells, or dry 
months 3,20 22 , although ref. 23 found the length of dry spells have 
increased over the conterminous US concurrently with observed 
increasing trends in precipitation intensity. 

Showing how a warming climate will affect dry- day frequency will 
help elucidate the processes involved in overall changes in annual 



precipitation and allow societies to better adapt to and mitigate the 
impacts of local hydrological changes. 

Results 

Changes in dry-day frequency. This study focuses on projected 
changes in dry-day frequency in a warming climate. We use daily 
precipitation from 28 climate models (all those with the required 
data that were available at the time of analysis) that participated in 
the CMIP5 project 24 . We take 1960-1989 as our base historical 
climatology period, and assess changes projected for the future 
(2060-2089). The historical simulations were forced with observed 
long-lived greenhouse gases (GHGs), solar radiation, and volcanic 
aerosols appropriate to the twentieth century, while the future 
simulations were forced with the Representative Concentration 
Pathway 8.5 (RCP 8.5) scenario. We also compare model 
climatology over the base period to observations from the daily 
version of the Global Precipitation Climatology Project data set 
[GPCP Daily; Ref. 25] during the restricted 1997-2012 period of 
GPCP Daily data availability. See the Methods section for further 
details. 

Partly due to coarse spatial resolution, global climate models pro- 
duce too many days with small precipitation totals 26 , sometimes 
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Figure 2 | CMIP5 multi-model ensemble average mean change in frequency of dry days (days/year) by 2060-2089, relative to the historical period 
1960-1989, using the RCP8.5 forcing scenario. Stippling indicates areas where at least 70% of the models agree on the sign of the change. Graph to the 
right: zonal means values. Map was produced using NCAR Command Language (NCL; Ref. 30). 



known as the "drizzle problem". Accordingly, we define days with 
precipitation < 1 mm as dry days for both models and observations 
(see Methods for further details). 

Fig. 1 shows the multi-model ensemble mean climatological 
(1960-89) yearly frequency of dry days in the historical simulations 
(panel a) and GPCP observations (panel b). The multi-model 
ensemble mean pattern is in good agreement with observations on 
broad spatial scales (spatial correlation of 0.91), capturing the max- 
ima in subtropical high-pressure regions and in the polar regions, 
and minima in the tropical convergence zones and subpolar regions. 
However, regional differences between the models and observations 
remain. For example, high numbers of dry days in the subtropical 
subsidence regions systematically extend farther west in the observa- 
tions than in the ensemble mean. Over most land areas, model dry- 
day frequencies are within 1 to 10 days of the observed pattern, 
except over tropical Africa, the Tibetan plateau, and subtropical 
South America. Over most ocean areas, the ensemble mean tends 
to underestimate the number of dry days, with errors in some loca- 
tions reaching —80 days/year. High frequencies (>300 days/yr) of 
dry days are observed over most of the arid subtropical regions, such 
as the Sahara Desert, Afghanistan, part of northwest India, Pakistan, 
the Southwest USA, parts of Peru, Chile, South Africa, Australia, and 
the polar regions. Models underestimate the dry-day frequency over 
these regions by 1-5 days. Recomputing Fig. 1 over the common 
period of the models and observations, 1997-2005, makes almost 
no difference (online supplementary material Fig. SI). 

The ensemble-mean dry-day frequency changes (future minus 
historical) in response to anthropogenic forcing are shown in 
Fig. 2. Increases in dry days are nearly all located in the extratropical 
regions between 50°N and 50°S; the Amazon basin and western 
Indonesia are the only tropical regions that experience significantly 
increased numbers of dry days per year. A strong increase in dry-day 
frequencies is projected over land regions around the Mediterranean 
and Caribbean Seas, Central America, Mexico, Chile, and western 
Indonesia. Increases are also projected over the oceanic regions west 
of Mexico, the eastern South Pacific and North Atlantic, and south- 
west of Australia. Some of these regions are projected, in ensemble 
mean, to experience an additional 30 dry days per year by the second 
half of the 21 st Century. Reduced numbers of dry days are projected 
in the intertropical convergence zone and polar regions. The 



majority of models agree on the sign of dry-day frequency changes 
over the regions with strongest projected increases or decreases (see 
stippled regions in Fig. 2), however, individual models diverge con- 
siderably in the magnitude of change. 

The projected increases in numbers of dry days in the region 
around the Mediterranean Sea extend east to the Caspian Sea, and 
represent at least a 10% increase in dry-day frequency, reaching 30% 
over the northern Iberian Peninsula. An additional 15-20 dry days 
per year (5-8% increases) are projected to occur in South Africa, 
Mexico, and western Australia. The most frequent (not the most 
intense) precipitation events in these Mediterranean climate regimes 
are associated with the equatorward flank of midlatitude cyclones; 
these are the events that will become less frequent with a poleward 
shift in the midlatitude stormtrack 10 . On the other hand, the events 
associated with the most intense precipitation usually follow more 
equatorward tracks; these extremes will still occasionally occur and 
should yield more precipitation in a warmer world 1415 . Over the 
Northern polar regions the number of dry days declines by about 
40 days per year (12%), with a smaller decline over the Southern 
polar region of —10 days per year (5%). On the other hand, over 
the currently dry subtropical regions of the Sahara Desert, 
Afghanistan, and northwestern India the models show little change, 
which is consistent with poleward expansions of subtropical dry 
regions and poleward shifts of storm tracks and jet streams 10-1113 . 

Contributions of dry days to annual precipitation changes. Future 
increases or declines in annual precipitation are an aspect of climate 
change that will have a particularly strong effect on societies and 
ecosystems. The multi-model ensemble change in annual precipi- 
tation projected by the models used here is shown in Fig. 3a, and is 
consistent with previous results based on CMIP3 and CMIP5 
models 3-5 . 

The question we ask is what fraction of the annual change can be 
attributed to changes in the number of dry days, as opposed to 
changes in precipitation intensity on wet days? Local impacts could 
be rather different depending on which mechanism dominates. We 
estimate this partitioning with the linearized approach of ref. 27 (see 
Methods for details); results are shown in Fig. 4. For comparison, 
Fig. 3b shows the change in precipitation intensity on days when 
precipitation occurs. 
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Figure 3 | CMIP5 multi-model ensemble average percent change in (a) annual mean precipitation; (b) precipitation intensity during precipitating 
days. Values are computed over the period 2060-2089, relative to the historical period 1960-1989, using the RCP8.5 forcing scenarios. Stippling indicates 
regions where at least 70% of the models agree on the sign of the changes. Graphs to the right: zonal mean values. Maps were produced using 
NCAR Command Language (NCL; Ref. 30). 
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Figure 4 | CMIP5 multi-model ensemble average percentage contribution of changes in dry days to changes in total annual precipitation. Values are 
computed over the period 2060-2089, relative to the historical period 1960-1989, using the RCP8.5 forcing scenario. Negative values show where 
changes due to dry days and precipitation intensity are competing. Map was produced using NCAR Command Language (NCL; Ref. 30). 
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Figure 5 | CMIP5 multi-model ensemble average change in (a) interannual variability of precipitation scaled by historical mean annual 
precipitation; (b) frequency of precipitation days (%). Values are computed over the period 2060-2089, relative to the historical period 1960-1989, using 
the RCP8.5 forcing scenarios. Maps were produced using NCAR Command Language (NCL; Ref. 30). 



Annual precipitation changes are determined by the interplay 
between changing frequencies of dry days (Fig. 2) and precipitation 
intensity on wet days (Fig. 3b). An example helps make this clear. In 
the Mediterranean Sea region, annual precipitation is projected to 
drop by 20% (Fig. 3a). Why does this happen? Fig. 3b shows that the 
intensity of precipitation on wet days is projected to decrease in this 
area, but only by about 5% and not consistently across models. By 
contrast the number of dry days is projected to increase more than 20 
days/year, a result the models agree on (Fig. 2). The net result is that 
more than 80% of the projected annual precipitation change in the 
Mediterranean Sea region can be ascribed to increases in occurrence 
of dry days (Fig. 4). 

Increases and decreases in dry-day frequency account for more 
than half of the change in total-annual precipitation in regions lying 
between 40°N and 40°S (Fig. 4). In particular, Figs. 3 and 4 show that 



declines in annual precipitation over subtropical regions are caused 
primarily by increases in dry-day frequencies. By contrast, moisten- 
ing of the tropics and higher latitudes, and drying of the equatorward 
flank of extratropical storm tracks, are driven by changes in both 
frequency of dry days and precipitation intensity. In other words, 
model projections of increased aridity over subtropical regions are 
mostly due to increased numbers of dry days and are much less a 
result of declining average precipitation intensities. In the Amazon 
and Indonesian regions, the large increases in dry days projected to 
occur (Fig. 2) are the dominant factor accounting for annual precip- 
itation changes (Fig. 4); however the magnitude of the projected 
annual changes (Fig. 3a) is small, because the climatological numbers 
of dry days are low in those regions (Fig. 1). These increases in dry 
day frequency over the Amazon and Indonesian regions could be due 
to changes in ENSO teleconnections in a warmer climate 28 . 
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Figure 6 | CMIP5 multi-model ensemble average change in frequency (days/year) of days with precipitation in the (a) 5-25th (b) 50-75th (c) > 99th 
percentiles. Maps were produced using NCAR Command Language (NCL; Ref. 30). 
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In addition to their importance to changes in mean annual pre- 
cipitation, the changes in dry day frequency also affect interannual 
precipitation variability. In essence this happens through changes in 
sample size: more dry days means fewer precipitation events con- 
tributing to the annual total. Smaller annual sample of days with 
precipitation then leads to more variability about the annual mean. 
This is illustrated in Fig. 5, where the model-projected change in 
interannual precipitation variability (Fig. 5a, calculated as the stand- 
ard deviation of annual total precipitation normalized by the his- 
torical mean) is compared to the change in the frequency of 
precipitating days (Fig. 5b). The patterns differ in polar regions, yet 
globally, the spatial correlation between the two is —0.79, reflecting 
the close link between anthropogenic trends in daily precipitation 
frequency and interannual variance. As a sensitivity test, we addition- 
ally calculated the interannual variability as the difference between 
the 95 th and 5 lh , 90 th and 10 th , or 75* and 25 th percentile values of total 
annual precipitation instead of it's standard deviation and found the 
same result. The association between increasing numbers of dry days 
and greater interannual variability is strongest over Mediterranean- 
type climate regions where rainfall is infrequent, strongly seasonal, 
and already highly variable. In such regions an anthropogenically- 
driven increase in the numbers of dry days will increase the already 
large natural interannual precipitation variability. 

The changes in dry-day frequencies reported here are part of 
broader shifts within frequency distributions of days with a given 
precipitation rate [cf. Ref. 3, who examined annual mean accumula- 
tion in precipitation categories]. Fig. 6 shows the multi- model 
ensemble mean change in the frequency of days with precipitation 
in the 5-25 ,h , 50-75 lh , and >99th percentiles (calculated with respect 
to the historical period). Importantly, there is a projected broad scale 
tendency over many global regions towards a decline in the number 
of days/year with daily precipitation below the 75 lh percentile thresh- 
old (Fig. 6a,b). The polar regions and deep tropics are exceptions, 
where there is an increase in the frequency of days with light to 
moderate (5 th -75 lh percentiles) precipitation. By contrast, frequen- 
cies of days that experience very heavy precipitation (>99 th percent- 
ile) increase over most of the world, excepting most oceanic 
subtropical high pressure regions (Fig. 6c). In polar regions and 
the deep tropics, frequency of all precipitation categories increases. 
In sum, over much of the globe the numbers of dry days increase, and 
the intensity of extremely wet days increases, both at the expense of 
light to moderate precipitation days. 

Discussion 

Daily precipitation projections from 28 CMIP5 global climate mod- 
els forced by the RCP 8.5 emissions scenario are used to evaluate 
changes in dry day frequency and its effects on annual precipitation 
mean and variability during the period 2060-2089 relative to 1960- 
1989. In most subtropical semi-arid regions the number of dry days is 
projected to increase markedly (10-15%). These dry-day increases 
are particularly pronounced around the Mediterranean Sea, Central 
America and Mexico, the Amazon basin, Chile, and western 
Indonesia. A linearized analysis indicates that projected annual mean 
precipitation changes in those regions are dominated by an increase 
in the number of dry days rather than by changes in precipitation 
intensity. This anthropogenically-forced drying will have important 
hydrological, ecological and societal effects. The greater dryness will 
be aggravated by increased year-to-year variability, particularly in 
Mediterranean-type climate regimes where annual precipitation is 
already highly variable, since more dry days means that fewer wet 
days will contribute samples to the annual mean, resulting in greater 
interannual sampling variability. 

In some regions, such as the U.S. Southwest, the projected increase 
in dry days is more robust and consistent across models than the 
annual precipitation trend, since the latter reflects a competition 
between more dry days and greater precipitation intensity on wet 



days 27 . Such anthropogenic trends in the daily precipitation regime 
suggest increasing regional probabilities of drought as well as flood, 
but are obscured by temporal aggregation to multi-year averages. A 
better mechanistic understanding of dry-day occurrences and their 
geographic shifts will require regionally focused investigation, given 
the complex, spatially varied interplay of the influence of dry and 
wet extremes in shaping annual precipitation volumes and their 
variability. 

Methods 

Data. The daily model precipitation data were obtained from 28 coupled general 
circulation models that participated in CMIP5 experiment 24 , available at http:// 
pcmdi9.llnl.gov/. A list of the models used can be found in the online supporting 
material. We used historical (20 th century) and Representative Concentration 
Pathway 8.5 (RCP 8.5) scenarios to represent present and future climatology, 
respectively. The historical simulations were forced with observed long-lived 
greenhouse gases (GHGs), solar radiation, and volcanic aerosols for the observational 
period, while the RCP 8.5 simulations were forced with an increase of the atmospheric 
concentration of C02 to 950 ppm by year 2100. 

The observed precipitation data set used is the GPCP Daily version 2 (Ref. 25), 
which covers the period 1997 to 2012 with a 1° X V spatial resolution. The GPCP 
Daily data is provided by the NASA Goddard Space Flight Center, Mesoscale 
Atmospheric Processes Laboratory, from their website at http://precip.gsfc.nasa.gov. 
For consistency, all modeled and observed precipitation was remapped to a common 
2° X 2" latitude-longitude grid by bilinear interpolation before analysis. 

Definition of dry day. Typically, the number of days with small precipitation 
accumulation is over-simulated in current GCMs, which have only relatively coarse 
spatial resolution 26,29 . Observations also do not report the precipitation value below a 
certain threshold; for example, NOAA's co-operative observing stations typically do 
not report values smaller than 0.1 mm/day. For consistency between models and 
observations we select a threshold of 1 mm/day, below which both models and 
observations re-gridded to the same grid are considered to have a dry day. The 1 mm/ 
day threshold is primarily dictated by the coarse resolution of the global models used, 
and is found to be the best overall choice after considering other possible values, 
including fixed thresholds of 0.1 mm/day and 0.5 mm/day as well as location- 
dependent thresholds consisting of the 5 th , 1 0 th , 1 5 th , 25 th , or 35 th percentile of the local 
precipitation on wet days. 

Effect of changes in dry days on annual precipitation. Following Ref. 27, a linear 
estimate of the effect on total annual precipitation of an additional dry day in a given 
month in the future is calculated by subtracting the amount of precipitation 
experienced on an average precipitating day in the same month in the historical 
period. The change in total annual precipitation that is due to changes in daily 
precipitation intensity on wet days is then taken as the total annual precipitation 
change minus the change that is due to dry day frequency. Although projected 
changes in intensity are not included directly in this approach, the largest projected 
contributions to drying are over subtropical regions where the smallest changes in 
intensity are projected. 

Graphic software. All maps were produced using NCAR Command Language 
(NCL). 
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